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Abstract 

The large N limit has been successfully applied to QCD, leading to qualitatively correct results 
even for iV = 3. In this work, we propose to treat the number TV = 3 of Standard Model generations 
as a large number. Specifically, we apply this idea to the neutrino anarchy scenario and study 
neutrino physics using Random Matrix Theory, finding new results in both areas. For neutrino 
physics, we obtain predictions for the masses and mixing angles as a function of the generation 
number TV. The Seesaw mechanism produces a hierarchy of order between the lightest and 

heaviest neutrino, and a ^13 mixing angle of order 1/N, in parametric agreement with experimental 
data when N goes to 3. For Random Matrix Theory, this motivates the introduction of a new 
type of ensemble of random matrices, the "Seesaw ensemble." Basic properties of such matrices 
are studied, including the eigenvalue density and the interpretation as a Coulomb gas system. 
Besides its mathematical interest, the Seesaw ensemble may be useful in random systems where 
two hierarchical scales exist. 
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1 Introduction 

The large limit is a very useful tool in various theoretical models, both at a qualitative and a 
quantitative level [1]. Historically, one of the most important examples came from 't Hooft's large 
analysis of QCD [2j. Although the number of colors Ac = 3 is not very large, many features of 
mesons and baryons calculated in a 1/Ac = 1/3 expansion are in good agreement with experiment \3\. 
One cannot help but notice that in the Standard Model (SM), the number "3" also appears as the 
number of fermion generations. Here, we propose to treat the number of families Nf = 3 in a large N 
approximation, and study the properties of fermion masses and mixings as a function of the expansion 
parameter 1/Nf = 1/3. 

In this paper, we make the first attempt to use the large A limit to understand some part of the 
flavor problem (i.e. the origin of masses and mixings of the SM fermions). More specifically, we will 
apply this idea to the SM neutrinos assuming the mass anarchy scenario of [U [5] , where the neutrino 
couplings are treated as random variables. In light of the recent measurement of ^13 from Daya Bay [6] 
and later confirmed by Reno [7] , the anarchy hypothesis is consistent with the measured value [H [9] . 
We will not consider the quark sector here (which requires additional structure), although we will 
comment briefly on possible large A applications in ^ 
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Extending the 3x3 random matrix for the three active neutrinos to an A'^ x random matrix 
(hereafter N = Nf denotes the number of generations), we have to analyze the statistical properties 
of the neutrino spectrum and mixing matrix. This is a classic problem in Random Matrix Theory 
(RMT); see e.g. [10\ [TT] for reviews. RMT has an incredible range of applications, from medical HIV 
studies \T2] to the string landscape [131 El [IS], and now we are adding neutrino physics to this list. 

RMT provides analytic results for the correlation functions of eigenvalues of large rank matrices. 
For example, the eigenvalues of Hermitian complex matrices obey the well-known Wigner semicircle 
distribution. Knowing the spectral density for the neutrino masses would already provide valuable 
information. It will allow us to calculate the expectation value of the ratio of the lightest and heaviest 
masses, which will in turn determine how hierarchical the mass spectrum is. Although experimentalists 
have measured all three mixing angles of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, it 
is still not clear whether the neutrino mass spectrum is normal or inverted, and hierarchical or quasi- 
degenerate. Our analysis for general A^ will reveal that the statistical properties of the spectrum 
depend sensitively on the underlying mass mechanism. 

There exist different possibilities for the neutrino mass matrices. Our main focus will be on the 
Seesaw mechanism, which provides a dynamical explanation for the smallness of neutrino masses. We 
will demonstrate that the Seesaw mechanism is not only required to explain the hierarchy between the 
neutrino mass scale and the weak scale, but also preferred by the large A^ behavior of our RMT analysis 
and the current experimental results on the neutrino mass-squared differences. Besides providing a 
useful tool to analyze the spectrum, the connection with RMT will also lead us to define a new type 
of random matrix ensemble, as we describe in more detail below. 

We have attempted to make our paper accessible to both the neutrino physics and RMT com- 
munities. We begin in ^ with a short review of the results in RMT that we will need, including 
the Gaussian unitary ensemble in §2.11 and the Coulomb gas picture in §2.21 Next, we apply RMT 
techniques to the simplest case of complex Major ana neutrinos in ^ Our main results are presented 
in ^ where we analyze the large A^ spectrum of random Seesaw neutrinos, and also define and study 
the new "Seesaw ensemble." In ^ we apply these results to the phenomenologically relevant A^ = 3 
case and compare with experimental data. Our conclusions and future directions are summarized in 

m 

2 Review of random matrix theory 

Let us begin by reviewing some basic results from random matrix theory that will be needed in the 
following sections. For a pedagogical exposition, we refer the reader to [TT]. 
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2.1 Gaussian unitary ensemble 

Random matrix theory is the study of the statistical properties of eigenvalues of very large matrices. 
We will first consider N x N Hermitian complex matrices M = A + , where each element of A is 
an independently distributed random variable. The space of such matrices is known as the Gaussian 
unitary ensemble (GUE). 

Writing the probability distribution as P(AI)dM , the defining properties of the GUE are that: a) 
it is invariant under unitary transformations, P{M)dM = P{M')dM' , with M' = U'^ AdU; and b) the 
matrix elements are statistically independent, namely P{M) = Y\i<^j Pij (Mij) . From this, it follows 
that [lO] 

P{M)dM oc dMe-f ^'■l*^^*^] , (2.1) 

where M has been linearly redefined to set the center at zero and to fix the coefficient in the expo- 
nential. The matrix M is diagonalized by 

M = UMJ^ . (2.2) 

Here, A = diag(mi, . . . , m^r) with rrii as the eigenvalues of the matrix with the convention rrii < rrij for 
i < j and U a special unitary matrix with N{N — l)/2 independent parameters. Changing variables 
to A and U then gives 

P{M)dM = Cat dC/ JJ drui JJ(mi - rrijf e'^ ^« . (2.3) 

i i<j 

Here, dU is the invariant Haar measure for U{N) and Cn is a normalization constant. 

We see that at large the distribution of angles and phases contained in U becomes flat in the 
natural group theory coordinates. This turns out to have important consequences for the neutrino 
mixing angles, analyzed by ^ i5j. On the other hand, the eigenvalues have nontrivial statistical 
properties. These are characterized by the k-point correlation functions 

A^! f 

Rk{mi, . . . ,mk) = j-z p- dmk+i...dmNP{mi,...,mN). (2.4) 



{N - k)\ 

Of particular interest for us will be the one-point function 



{p{m)) = l^^^5{m — rrii)^ = N J dm2 ■ ■ ■ drriN P{m, m2, ■ ■ ■ , . (2-5) 

This gives the eigenvalue density, and its inverse l/p{m) determines the mean level spacing around 
m. 
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The GUE can be solved exactly in the limit iV — ?■ oo. This is done by recognizing that the distri- 
bution in Eq. ()2.3p can be written as a product of orthogonal Hermite polynomials. The correlators 
are then of the form 

Rk{xi, . . . ,Xk) = det[KN{xi,Xj)]ij=i^,,,^k , (2.6) 
where the kernel K]^ is defined as 

Af-l 

Kn{x, y) = e-(^'+^')/2 Hn{x)Ka{y) • (2.7) 

n=0 

Here we have introduced the rescaled variables 

Xk = ]l—mk. (2.8) 

This result implies that at large the eigenvalue density is given by the Wigner semicircle distri- 
bution, 

'iV2iV - x2 forx2<2iV, 



{p{x)) 



(2.9) 

for x^ >2N . 



Changing back to the original eigenvalues m yields p{m) = 4 — for rrP < 4 and zero otherwise. 
Thus, the spacing of eigenvalues near the origin is of order Of course, the range of the semicircle 

(4 in this case) is conventional, being related to the normalization of the exponential in Eq. ()2.3p . 

2.2 The Coulomb gas picture and generalizations 

We will now review the physical interpretation of the GUE as a Coulomb gas, introduced by Dyson [ 16^ 
I17j . This picture will be useful when we study the "Seesaw ensemble" in ^ 

The basic idea is that Eq. (12. 3p can be interpreted as the thermal partition function for N charged 
particles with positions (xi, . . . ,xn) moving on a fixed line in two dimensions, under the influence of 
a harmonic potential and electrostatic repulsion. In more detail, the Hamiltonian for this system is 

i i<j 

Then the partition function 

Zn{(3) = Zo I JJdxie-^^ (2.11) 

i 

reproduces the GUE distribution if the inverse temperature is fixed at /3 = 1. [These expressions are 
given in terms of the rescaled variables in Eq. ([27 
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Physically, the properties of the eigenvalue distribution appear from a competition between the 
harmonic potential and the electrostatic repulsion. If the spacing between charged particles decreases 
with N, at large N the system admits a continuum approximation 



H = J dx p{x)x'^ — J dxdy p{x) p{y) log \x — y\ , (2-12) 



where p{x) is the level density introduced in Eq. (j2.5p l^ The density is determined by extremizing H, 
subject to the constraints p{x) > and J dxp{x) = N . The stationary condition leads to the integral 
equation 

— x^ + j /9(y) log |a; — ?/| = const . (2-13) 

The solution to this equation is given by the semicircle distribution Eq. ()2.9p . 

This picture suggests a way of understanding non-Gaussian ensembles, with probability distribu- 
tion 

P(M)oce-''^[^(*'^'^'^')]. (2.14) 
Such an ensemble is equivalent to a system of charged particles with a Hamiltonian 

H = Y, V{x.i) - log{xi - Xjf . (2.15) 

i i<j 

The particles interact with each other via the logarithmic Coulomb repulsion, and moreover each of 
them is subjected to a nontrivial confining potential V{x). The level density p{x) of the ensemble 
can be obtained, as before, by going to a continuum limit at large N and extremizing the analog of 
Eq. (|2.12p with replaced by V{x). The result is an integral equation for p(x), 

— V{x) + J dy p{y) log \x — y\ = const . (2-16) 

This relation will be applied in 21 where knowledge of the mass distribution generated by the Seesaw 
mechanism will be used to construct a suitable confining potential V{x). 

3 Majorana neutrinos and 1/N hierarchies 

It is helpful to begin our analysis with the simplest case of neutrinos with random Majorana masses. 
Results for Dirac neutrinos are qualitatively similar. We consider neutrinos with Majorana 
masses 

Cm = -lMi^iyfCiyj + h.c., (3.1) 



^The factor of 2 from log (a; — y)^ was canceled against a factor of 1/2 from the requirement i < j in Eq. (|2.10[) . 
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where C is the charge conjugation operator. My is an N x N complex symmetric matrix. For the 
SM, = 3 and Cm comes from a dimension- five operator {HL)"^ . As in [11 [5], My is taken to be 
random. Instead of speciahzing from the start to = 3, we will determine the properties of neutrino 
masses at large and then take iV — t- 3. At large N the lightest neutrino will be found to have a 
mass suppressed by The results from this section will also help to develop intuition and tools to 

study in ^the more interesting and nontrivial random Seesaw mechanism case. 

We first analyze the statistical properties of random complex symmetric Major ana mass matrices. 
In order to find the probability distribution of eigenvalues, we perform a singular value decomposition 

My = UAU^, A = diag(mi,...,m7v) , UeU{N). (3.2) 

The measure becomes [5] 

dMy = dU JJmfcdmfc \ mf — rnj\ . (3-3) 

k i<j 

As reviewed in fj2l statistical independence of the elements {My)ij and invariance under unitary trans- 
formations implies that the probability distribution is Gaussian. Therefore, the complex Majorana 
neutrino masses are distributed according to 

P{My)dMy = CNdUW_mkdmk JWrnf - m]\ b-tE.'"? , (3.4) 

k i<j 

which is known as the Altland-Zirnbauer CI ensemble The eigenvalues occur in pairs =bmj, which 
will be important shortly. (This is also why the exponential contains an extra factor of 1/2 as compared 
to the GUE.) 

The statistical distribution of neutrino masses is equivalent to a Coulomb gas system with Hamil- 
tonian 

H = - '^xj - ^\og \xi -x]\-Y^ log \xi\ , (3.5) 

i i<j i 

where The last term encodes an electrostatic repulsion between the images ixj. This 

term dominates near the origin, forcing the density of eigenvalues p{x) to vanish as x — )■ 0. For 
X ^ A^^^/^, the density is a linear function, while for larger x it is well approximated by the semicircle 
distribution |18] . In more detail, the approximate distribution of neutrino Majorana masses is 



p{m) ~ < 



A^ \m\ for \m\ ^ , 

(3.6) 

N\/4: - for ■^\m\ <2 , 
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Figure 1: Qualitative spectrum for the complex symmetric case with a large A^. The location of the 
lightest, median and heaviest eigenvalues are shown from bottom to top. 



where we have dropped numerical constants and the range of the distribution (here set to \m\ = 2) is 
arbitrary. Exact formulas for the eigenvalue density and higher correlation functions can be obtained 
in terms of Laguerre polynomials [19j. 

Let us now explore the phenomenological consequences of these results. In the = 3 case we 
are interested in the spectrum of masses of the lightest, intermediate and heaviest neutrino. This is 
generalized to > 3 by looking at the expectation values of the lightest eigenvalue mi, the median 
m.(7v+i)/2 (or depending on whether is even or odd), and the heaviest rriN- Given the eigenvalue 
density p{m), these expectation values are obtained from 

k= dm p{m) . (3.7) 
JO 

For eigenvalue distributions that are approximately constant near the origin (such as the semicircle 
distribution of the GUE), mi/mjv ~ l//9(0), which is the spacing of levels at small m. For complex 
Majorana neutrinos this estimate does not apply directly, because the distribution Eq. (j3.6p vanishes 
at the origin. However, the more precise Eq. (|3.7|) shows that the level spacing is parametrically of 
the same order as that of the semicircle law. Similarly, the median is found to be ~ 1/2. 
In summary, we have 

TTT-TV A^ m-AT 2 

We have also checked these results numerically, finding the same type of behavior. The neutrino 
mixing angles and phases are encoded in the unitary matrix U of Eq. p.2p . At large A^ these variables 
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are distributed according to the Haar measure. For — t- 3 there are small deviations depending on 
whether the scan over neutrino matrices is done over a hypercube or a hypersphere tr(M^M) < 1. 
The phenomenological predictions of this measure were studied in detail in [H [5] . 

This type of parametric spectrum is not the one suggested by experimental data (reviewed in at 
least in the limit where large results give a good approximation to the = 3 case. For this reason, 
we next turn to analyze the statistical properties of masses generated by the Seesaw mechanism. 

4 New hierarchies and ensemble from Seesaw neutrinos 

We now focus on the more interesting and nontrivial Seesaw models, which have been used to explain 
the smallness of the neutrino masses with respect to the electroweak scale. We will find that starting 
from random masses for the light active and heavy singlet neutrinos, the Seesaw mechanism leads 
to hierarchies of masses that are parametrically different from those obtained in ^ Furthermore, 
this mechanism will motivate the definition of a new type of random matrix ensemble, the "Seesaw 
ensemble" , whose interesting mathematical properties we will begin to explore here. 

4.1 Random Seesaw neutrinos 

Let us consider the scenario with N active (left-handed) neutrinos, and N singlet (right-handed) 
neutrinos which have heavy Majorana masses. The case with different numbers and Nji of left- 
and right-handed neutrinos also appears to be interesting, offering the possibility of taking the large 
Npt limit while keeping Nl fixed. We will comment more on this possibility in ^ The neutrino mass 
matrix is a 2N x 2N matrix of the form 



where Mpt is the symmetric complex Majorana mass matrix for the right-handed neutrinos, and 
Md is the complex Dirac mass matrix between the left-handed and right-handed neutrinos. For 
||M^|| ^ II-^dII the heavy SM singlet neutrinos can be integrated out to obtain the N x N matrix 
for the active neutrinos. 



This is the matrix that will be the object of our study. 

Before proceeding, it is necessary to point out that the overall scales for Mij and M/j will not be 
addressed in the paper. We only assume that ||M£)|| <C ||M/j||, so that the approximation Eq. (14.20 is 
valid and the Seesaw mechanism can be used to explain the smallness of the neutrino mass scale. 




(4.1) 



= Md M^^ Ml . 



(4.2) 
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We take M^) and M/j as random matrices and study the statistical properties of the eigenvalues 
of the matrix in Eq. (14. 2p . Starting from these ensembles and using Eq. (|4.2[) . it should be possible 
to determine the spectrum and correlation functions for the eigenvalues of My. However, such a 
first-principle derivation turns out to be very involved, requiring the diagonalization of and the 
Jacobian matrix to rewrite the distributions for M/j or Mpi in terms of the eigenvalues of M^. Instead, 
we will consider an alternative strategy based on an ansatz for the eigenvalue distribution p{mu) for 
My, which will be verified numerically. Then we will apply the Coulomb gas picture to write down an 
action for the eigenvalues m^, that can be used to derive the correlation functions of the system. 

The first step is to diagonalize M^ and Mr separately, rewriting My as 



My 



(mi 



V 



7712 



\ (Mr' 



niN ) 



Mo 



V 



M-nI 



(mi 



\ 



1712 



\ 



niN I 



(4.3) 

Here, U^ceX is a unitary matrix that comes from the relative diagonalization of Md and M^ and 
(Mi) are the eigenvalues of Md (M/j). The eigenvalues of My will be denoted by rriy. In order to 
develop intuition on the spectrum of My, let us approximate the eigenvalues of M^ and Mn by 



mi 



N 



■ m 



and 



Mi^—M. 
N 



(4.4) 



This is a good approximation away from the edges of the semicircle distributions. 

We will find that the distribution of My can be described in terms of an action Eq. (|2.15p . with a 
confining potential V{my) and a logarithmic Coulomb repulsion. Let us now consider the effects from 
eigenvalue repulsion. 



repul 



El 

i<j 



log {mi -ml) . 



(4.5) 



The matrix C/j-ei rotates the relative orientation of Mj and mj for the final eigenvalues. Although there 
are many possible relative orientations, we choose two points for illustration purpose. The first case 
is to have 



Case A : '^La = 

which has C/rei = ^N- The second example is 

Case B : "t-^b 



m^ 
M 



m? i 

'm N 



2 

mf 



m 



Mn+ 



l-i 



M N{N + 1 



(4.6) 



(4.7) 
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Figure 2: Qualitative spectrum for the Seesaw case with a large A^. The location of the lightest, 
median and heaviest eigenvalues are shown from bottom to top. 



with [/rei as a permutation matrix that interchanges ioA^ + l — zin Mr. In other words, this case 
corresponds to having the largest suppression for the smallest My eigenvalue. 

Comparing the repulsion Hamiltonian for these two choices of orientations, we have 



^repul,A - ^repul,B = " X] '^"^^ " " 



•2 \ 2 / -2 \ 2' 



N+l-jJ \N+1 



>0, (4.8) 



which shows that the case B has a smaller energy from contribution of the Coulomb repulsion. Con- 
centrating then on the eigenvalue distribution of the case B (we drop the label "i?" from now on), we 
have, for the smallest eigenvalue and the median. 



1 1 N/2 . 

ml 1 niy 1 



(4.9) 



This is shown qualitatively in the spectrum in Fig. [2j Based on numerical results, we will now argue 
that the case B gives a very good approximation to the full eigenvalue distribution. 

The spectrum of My is studied numerically as follows. Treating all elements of Md and Mr as 
complex numbers, we choose random numbers for both the real and imaginary parts of the elements 
of M£) within (—1,1) and (— lO^'', lO^'^) for M/j, which incorporates the two hierarchic scales in the 
Seesaw mechanism. We calculate M^, and diagonalize it. This procedure is repeated several thousand 
times to calculate the eigenvalue distribution. Let us now focus on the smallest eigenvalue and the 
median, and compare with the predictions Eq. ()4.9p . The full distribution will be studied below. 



10 



The numerical results for the smallest eigenvalue as a function of A'^ are shown in the blue circles 
of the left panel of Fig. [3l together with a fitted distribution using a function proportional to N~'^. 
One can see that the is indeed a good description of the large N behavior and also of the = 3 
case. The right panel of Fig. [3] shows the ratio of the median eigenvalue over the largest eigenvalue. 
(For simplicity, only odd N numbers are shown in this plot.) The numerical results are well fitted by a 
N^^ function. Therefore, we find a good agreement between Eq. ()4.9p and the numerical eigenvalues; 
this provides nontrivial evidence that the distribution Eq. (j4.7p is a good approximation. Further 
properties of the eigenvalue distribution will be analyzed shortly. 




Figure 3: Left panel: the blue circle points are the ratios of the smallest eigenvalue over the largest 
eigenvalues as a function of the matrix rank for the Seesaw case. The red continuous line is the fitted 
result using a function proportional to N^^. Right panel: the same as the left but for the median 
eigenvalue over the largest eigenvalue. The fitted function is chosen to be N~^. 



Since we are doing statistical analysis for the expectation values of the eigenvalue ratios, it also 
important to compute the standard deviation of these observables as a function of A^. Defining 




(4.10) 



the large dependence of (Tmi/nij turns out to be the same as that of mi/rrij. This is illustrated in 
Fig. m which shows the numerical ratios of the variance over the expectation value. This implies that 
our large A^ analysis can not capture detailed order one coefficients in front of the A^ dependence. 

Comparing Fig. [1] and Fig. [2] shows that the Seesaw mechanism generates a larger hierarchical 
spectrum than the one obtained from random Majorana masses. For a fixed heaviest eigenvalue, 
random Majorana neutrinos have a lightest eigenvalue suppressed by 1/A^, while in the Seesaw case 
the suppression is by l/N"^. The behavior of the median is also parametrically different as a function of 
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Figure 4: The ratios of the variance over the expectation value as a function of A^. 



A^. These statistical properties are then very sensitive to the mass mechanism, providing predictions 
that distinguish the two scenarios. In ^ we will compare those two spectra to the experimentally 
observed values, and conclude that the Seesaw scenario is preferred. 

4.2 The Seesaw Ensemble 

So far we have found that the Seesaw mechanism provides an interesting way of generating additional 
hierarchies out of randomness. Starting from a small input parameter of it produces a parametri- 
cally smaller l/N"^ suppression. In terms of the original masses in Mr and Mjj, with ||Mjj|| ^> ||M£)||, 
this is the statement that the expectation values are dominated by the relative orientation Eq. ()4.7p . 
This effect may also have applications in other systems. For this reason, we now study in more detail 
the "Seesaw ensemble" by itself, independently of motivations from neutrino physics. As far as we 
know, this provides a new type of ensemble beyond the ones studied in the literature (see e.g. [IT] for 
a recent review). 

We define the Seesaw ensemble as the set of N x N matrices of the form 

M^ = MdM^^mI, (4.11) 

where Mq and Mn are complex random N x N matrices, and Mr is symmetric. Recalling the 
discussion in iSl Mr belongs to the AZ-CI ensemble^ On the other hand, the probability distribution 
for the Dirac matrices M^ is given by [5] 

P{MD)dMD oc JJmfedmfc fj(mf - m]f e-^^-"^^- . (4.12) 

k i<j 

^In particular, the eigenvalue density vanishes at the origin. 



12 



This is in the same class as the unitary ensemble that appears in QCD, with no flavors and a nonchiral 
Dirac spectrum. 

Ideally we would want to derive the probability distribution for My starting from those of M^, 
and Mr. This is, however, a very complicated task. Instead, we now explore the consequences of the 
ansatz Eq. ()4.7p for the spectrum of eigenvalues of My, which we found to be in good agreement with 
numerical results. Our analysis of the Seesaw ensemble will rely heavily on the generalized Coulomb 
gas picture described in ^2.2[ 

From Eqs. (|3.7p and (j4.7p . it is not difficult to derive the continuous spectral density for the Seesaw 
ensemble, 



^ . ^ _ 1 xN'^ + 2N + 2-N Vx^N^ + 4xN + 4x 

2 iV + 1 \/x^N^ + AxN + 4x 

where the normalization is such that dxp{x) = N. For very small and large values of x, the density 
behaves as 

_/Y2 + 1 

x^O: PSees^y^ix) 2VN + ' X ^ oo : Psees^^ix) . (4.14) 

Both limits reveal interesting properties of the ensemble. Note that the spectral density is divergent 
when a; — )• 0, but p{x) is still integrable. Furthermore, p{x) does not have compact support -there 
is no upper bound on eigenvalues for the Seesaw ensemble. This comes from the inverse of Mr in 
Eq. ()4.1ip . so that it is always possible to have eigenvalues much larger than the light scale m? /M. 
These properties are in sharp contrast with the GUE or AZ-CI ensembles that we discussed before, 
for which p{x) has compact support and is finite as x — )• 0. 

In order to verify our analysis so far, let us compare the spectral density Eq. (I4.13P with numerical 
results. This is done in the left panel of Fig. [5l which shows both the numerically averaged density 
of eigenvalues for N = 200 and the analytic prediction. One can see a good agreement between those 
two distributions. 

Let us now develop the Coulomb gas model for the Seesaw ensemble. As we reviewed in f|2l this 
representation is given by a system of particles moving along a line in two dimensions, subject to a 
Coulomb repulsion and a confining potential V{x). We already discussed the effects from repulsion, 
and now we will determine V{x). Plugging the eigenvalue density into Eq. p.l6p obtains 



Seesaw {x) 



2{N + 1) 



_/Y + l\ / / N'^x 

Nx log ( — + 2\/ N'^x'^ + ANx + 4x tanh"^ 



N'^x I V ' ' 1 Y + 4iV + 4 



(4.15) 
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Figure 5: Left panel: the normalized spectral density as a function of eigenvalues for the Seesaw 
ensemble. The points for the analytic results are calculated by integrating p{x) for the bin size. Right 
panel: the red (solid) and blue (solid) lines are the effective confinement potentials for N = 100 and 
= 200, respectively, in the Coulomb gas picture. The black (dashed) line is the harmonic potential 
for the CUE. 



Mapping back to the matrix eigenvalues, this means that the probability distribution for Seesaw 
matrices is 

P{M^)dM^ (xYldnikYl \m\ - m|| e-^^^*"?) . (4.16) 

k i<j 

In the limits of small and large values of x, we have 



■ (4.17) 



Thus, the confining potential V'g^gj,3^^(x) is continuous at = 0, but its first derivative diverges there. 
This is much steeper than the harmonicpotential for the GUE, explaining why the eigenvalue density 
is parametrically larger near the origin^ The Coulomb gas picture hence provides a physical way of 
understanding the steeply decreasing behavior of p{x) as x increases. 

Although we will not study higher correlators in detail, let us mention that given V{x), the 
correlation functions for the Seesaw ensemble can be obtained as functions of polynomials that are 
orthogonal with respect to e~^^^\ They are of the form Eq. (12. 6p in terms of these new orthogonal 
polynomials. 



^See for instance [20] for different potentials in other q-random matrix ensembles. 
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5 Application to neutrino physics 



Having analyzed the statistical properties of neutrinos for general N, let us now focus more specifically 
on neutrino physics, comparing our results (extrapolated to = 3) with the experimental data. 

From all neutrino experimental results including solar, atmospheric, reactor and accelerator exper- 
iments, the global analysis [21j obtains the following numbers for the neutrino mass squared differences 
and mixing angles, defined in the PDG |22j : 

Amii = (7.50 ± 0.185) x 10"^ eV^ , (5.1) 

Ami m = (2.47l°:r7) X 10"^ eV^ , Amh [I] = (-2.431°;°^^) ^ ^^-3 ^ (5 2) 

sin^ = 0.30 ± 0.013 , sin^ 623 = 0.41+[]:°^7 © 0.59t[!:[]^^ , (5.3) 

sin^ 6*13 = 0.023 ± 0.0023 , 6cp = 300° +i|*^8 . (5.4) 

Here, Am?j- = m? — m^; "N" means the normal ordering for the neutrino spectrum such that Am^^ <^ 
(Am|2 ~ Am|i > 0); "F means the inverted ordering for the neutrino spectrum such that Am^i ^ 
— (Amg-j^ ~ Am'^2 < 0). We take their fitted results for treating the reactor fluxes as free parameters 
and including the results for the reactor experiments with L < 100 m. As emphasized in [2T|, although 
a zero Dirac CP violation phase 6cp is referred, the significance is less than 1.5o" for the normal ordering 
and 1.75(7 for the inverted ordering. The absolute scale of neutrino masses is currently unknown and 
bounded from above as '^rrii, < 0.28 eV at the 95% confidence level [23], assuming a flat ACDM 
cosmology. 

We should emphasize that our analysis does not explain the overall scale of neutrino masses, but 
only ratios of masses. Let us concentrate flrst on the ratio Am'^i/Am'^^ ~ 0.03. For complex Majorana 
masses without the Seesaw mechanism (see Fig. [T|), the prediction is Amgi/Amgj^ = 0(1) at large 
A^, which does not provide a good fit to the experimental value. On the other hand, the prediction 
for Seesaw masses (see Fig. E]) is Aml^/Aml^ = 0{l/N^) = 0(0.1) for = 3, much closer to the 
experimental result. (Recall that these predictions have order one uncertainties.) Therefore, assuming 
random masses, the Seesaw mechanism is favored by data. 

Focusing then on the Seesaw mechanism, let us make the following observations based on the 
spectrum in Fig. [2j First, the neutrino spectrum is a normal hierarchical spectrum. The general 
neutrino mass scale is determined by the heaviest mass, ~ y^'Am^ ~ 0.05 eV, which could be tested 
from the Cosmic Microwave Background by the Planck satellite or the KATRIN experiment [23] in the 
future. The prediction for the effective Majorana mass | (m) | , which may be measured from neutrinoless 
double beta decay experiments, is in the region of (2.3 x 10^^, 5.0 x 10^^) eV for the normal hierarchical 
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Figure 6: Left panel: the absolute value of the off-diagonal element Uuy as a function of N. Right 
panel: the absolute value of the Jarlskog invariant multiplying N'^, which denotes the size of the 
physical phases in the unitary matrix, as a function of N. 



spectrum [22]. Secondly, the lightest neutrino mass is predicted to be ms/N^ = ms/S^ w 0.002 eV, 
which is unlikely to be tested in the near future unless one can come up a new experimental method 
to measure the individual neutrino mass more precisely. 

Finally, let us discuss the basic properties of the mixing matrix for this case. As mentioned in [5j, 
the basis-independence assumption predicts that the PMNS matrix follows the Haar measure of the 
Lie group. While in principle the corresponding measure for the Seesaw ensemble can be obtained 
starting from the Haar measure distributions for M/j and M^, for our purpose it is enough to evaluate 
numerically the expectation values of the elements of the matrix U that diagonalizes mIm^. The left 
panel of Fig. [6] shows the behavior of the absolute value of the off-diagonal entry Uin as a function 
of A^, which is proportional to \/^fN. This is just what we expect based on the normalization of 
eigenvectors. For the phases in the rotation matrix, we use the generalized Jarlskog invariants to study 
the rephasing invariant CP violation [25^ [26] . In the right panel of Fig. [6l we show the distribution 
of A^^ |Im[C/ii?723C/*3[/2i]| as a function of A^. The overall factor A^^ cancels the A^ dependence of 
the absolute values of the elements. We can see from this panel that the phase part of this Jarlskog 
invariant is order 1 and independent of A^. 

Translating these results into mixing angles and phases, our prediction for sin^is is order 1/A^ = 
1/3, not far from the measured values sin^ia ~ 0.15. For the CP violation phases, we anticipate all 
the phases in the leptonic mixing matrix in the weak charged current interactions to be order 1. The 
predictions of our analysis are summarized in Table [1] 
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Spectrum 


Normal hierarchical 


Lightest neutrino mass 


v/Am§2 xO[{N = 3)-3] « 0.002 eV 


sin 6*13 


O [{N = 3)-i] w 0.3 


Sep 


Oil) 



Table 1: The large predictions for the neutrino properties for the three active neutrinos. 

6 Conclusions and future directions 

In this work we proposed to treat the number = 3 of Standard Model generations as a large number, 
studying the consequences of this idea for neutrinos. Assuming anarchic neutrino masses, this becomes 
a problem in Random Matrix Theory. We analyzed the neutrino spectrum and mixing matrices 
analytically and numerically, and compared with experimental results. The statistical properties of 
these observables were found to be highly sensitive on the underlying mechanism for neutrino masses. 
In particular, the Seesaw mechanism leads to a hierarchy for the lightest neutrino, and has a 

spectrum that is qualitatively compatible with experimental data. 

Already at this level, the connections between neutrino physics and RMT appear to be interesting. 
Based on the Seesaw mechanism, we defined a new type of random matrix ensemble, the "Seesaw 
ensemble," which, as far as we know, has not been analyzed before in the literature. The properties 
of this ensemble are qualitatively different from those of Gaussian distributed matrices. We proposed 
an ansatz for the spectrum that was in good agreement with the numerical results. In the Coulomb 
gas picture, the confining potential for the Seesaw matrices is much steeper than in the Gaussian case, 
with a diverging first derivative at the origin. 

It will be interesting to determine the correlation functions of the Seesaw ensemble starting from 
the distributions of the underlying matrices and Mpt in Eq. ()4.1ip . We hope that the analytic 
and numerical results reported here can shed light on this problem. There are also various possible 
extensions. First, one could consider directly the ensemble of matrices of the form Eq. (j4.ip . without 
restricting from the start to the limit ||-/Vf/j|| ^ ||M£)||. Furthermore, there is also a natural two- 
parameter generalization of this ensemble, where Mr is a square Nfj^ x Nr matrix, but Md is an 
Nl X Nr matrix. Physically, this corresponds to having Nl active and Nr singlet neutrinos [27]. One 
could then study the statistical properties at fixed A''^ and large Nr. More generally, the existence 
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of the new parameter Nr/Nl can lead to different large A'" behaviors or multiscale effects. This is 
reminiscent of QCD with A'^^ fundamental fermions, where taking Nf/N^ ~ 1 led to the Veneziano or 
topological limit [28] . 

The large limit can also be applied to neutrino physics with one or more sterile neutrinos 
or the SM quarks and charged leptons. The dramatically hierarchical pattern of quark masses and 
angles implies that the couplings in the Lagrangian cannot be randomly distributed -more structure 
is required to explain the flavor problem. Nevertheless, extra-dimension models such as [29] still allow 
for random coefficients, and it would be interesting to understand how the dynamics in the quark 
sector is modified by large N statistical effects. 

Finally, given the incredible range of applications of RMT, the Seesaw ensemble may occur in other 
types of systems. One would need some pattern of the form Eq. (j4.1ip . arising from the interplay of 
two hierarchical scales. 
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